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We study the energetics and dispersion of anomalous dimers that are induced by the Pauli blocking 
effect in a quantum Fermi gas of majority atoms near interspecies resonances. Unlike in vacuum, we 
find that both the sign and magnitude of the dimer masses are tunable via Feshbach resonances. We 
also investigate the effects of particle-hole fluctuations on the dispersion of dimers and demonstrate 
that the particle-hole fluctuations near a Fermi surface (with Fermi momentum hkp) generally 
reduce the effective two-body interactions and the binding energy of dimers. Furthermore, in the 
limit of light minority atoms the particle-hole fluctuations disfavor the formation of dimers with a 
total momentum hkF, because near hkF the modes where the dominating particle- hole fluctuations 
appear are the softest. Our calculation suggests that near broad interspecies resonances when 
the minority-majority mass ratio ms/mF is smaller than a critical value (estimated to be 0.136), 
dimers in a finite-momentum channel are energetically favored over dimers in the zero-momentum 
channel. We apply our theory to quantum gases of 6 Li 40 K, 6 Li 87 Rb, 40 K 87 Rb and 6 Li 23 Na near 
broad interspecies resonances, and discuss the limitations of our calculations and implications. 



I. INTRODUCTION 

Pairing of fermions or electrons in quantum condensed 
matter systems for long has been one of the most remark- 
able phenomena discovered in many-body physics [U-Q- 
It is also the fundamental issue in the so-called color- 
superconductivity phenomenon in the core of a neutron 
star Q. Recent developments in the field of ultracold 
atoms, especially the observation of atomic Feshbach res- 
onances, have further led to tremendous new opportuni- 
ties of studying novel pairing correlations in quantum 
gases [5-8]. Because of tunability of interactions and 
availability of a rich variety of alkali isotopes, cold gases 
have been turned into a marvelous, promising platform 
to study diversified pairing phenomena with unique cor- 
relations. Pairing properties near Feshbach resonances 
are mainly determined by, apart from scattering lengths, 
asymmetries in chemical potentials of the different atoms 
involved. Such asymmetries and consequently a very 
rich class of pairing phenomena are known either due 
to the imbalance in populations ffjl-flij or tunable mass 
ratios [l5l - tl8j . or both. Moreover, differences in the quan- 
tum statistics of underlying atoms also contribute to the 
asymmetries in chemical potentials as in Fcrmi-Bose mix- 
tures fl9l - [26| . Previous theoretical studies on Fermi- Bose 
mixtures have been focused on the formation of molecu- 



lar Fermi surfaces near narrow resonances 27J or far away 
from resonances [28j , pair correlations 2£ | , and quantum 
and thermal depletion of condensates 301 ] . 



A very closely related fascinating issue that has been 
challenging to the cold atom community is whether near- 
resonance two-body scattering in a channel with nonzero 
total momentum dictates the instability in interacting 
quantum mixtures. A major challenge in the studies 
of pairing phenomena in various quantum mixtures near 



resonances, unlike in the BCS-BEC crossover regime of 
a Fermi gas, is that there are more competing pairing 
schemes and possibilities of having higher order corre- 
lations; energetically it is quite difficult to differentiate 
competing scenarios near resonances. One natural ap- 
proach is of course to perform a full-scale numerical sim- 
ulation to resolve the issue. An alternative is perhaps to 
study few-body physics in the presence of a finite-density 
quantum gas and to gain insight on many-body correla- 
tions by exploring implications of few-body physics. A 
few interesting attempts have already been made along 
this direction. For instance, the properties of a single 
minority atom submerged in a quantum Fermi gas have 
been thoroughly investigated as a diagnosis of many- 
body correlations near interspecies Feshbach resonances 
[3ll - [38l ]. For Fermi-Bose mixtures, anomalous dimers 
with tunable masses were emphasized and the leading 
effect of particle-hole fluctuations had been studied dia- 
grammatically ;39]. At the same time, for Fermi- Fermi 
mixtures various instability lines have been determined 
numerically by studying the dimer and trimer formation 
in a truncated Hilbert space and the role of anomalous 
dimers as well as a universal trimer was explored [l8j . 
And very recently, Efimov states in a Fermi gas were 
studied in a few limiting cases and the spectrum flow has 
been obtained in the static- Fermi-sea approximation [40| . 
Logically speaking, these studies should form potential 
building blocks for constructing many-body states; they 
can also serve as a starting point for more systematic 
studies on the interplay between few- and many-body 
physics in cold gases near resonances. 

In this paper, we take a further step along this direc- 
tion hoping that our efforts to understand dressed bound 
states can shed more light on many-body pairing phe- 
nomena in mixtures. Particularly, we consider the limit 



2 



of a single atom or a minority atom in resonance with 
majority Fermi atoms which form a Fermi sea. This mi- 
nority atom can be either a fermion or boson although 
most of our discussions are in the context of a Bose atom 
resonating with a Fermi sea. We obtain the dispersion of 
dimers and investigate the Pauli blocking effect, the effect 
of fluctuating particle-hole pairs near the Fermi surface 
on the dimer energetics. These results should be appli- 
cable to quantum mixtures with extremely imbalanced 
population. We further examine the dimer dispersion 
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near broad interspecies resonances for 6 Li- 87 Rb 
6 Li- 23 Na [Hill, 40 K- 87 Rb HillH and 6 Li- 40 K 
that have been available in laboratories. 

Our analysis on the dispersion of dimers further sug- 
gests that when the minority atom is very light, minority- 
majority atoms near broad resonances would prefer to 
form dimers in finite-momentum channels that are ener- 
getically favorable. Although quantum mixtures of 6 Li- 
40 K, 6 Li- 87 Rb, 40 K- 87 Rb and 6 Li- 23 Na so far studied in 
laboratories are not in this particular limit, we antici- 
pate such a limit can be explored in future generations of 
quantum mixtures. More importantly, when combining 
the optical lattices and Feshbach resonances 0, EtJ , one 
can achieve resonances with continuously tunable widths 
and locations (48[. By properly choosing laser intensities, 
or by using isotope selective optical lattices (47[, the ra- 
tio of band masses of Fermi-Bose atoms can be further 
continuously tuned over a very wide range. This leads 
to potential opportunities of studying the limit of light 
Bose atom. 

The rest of the paper is organized as follows. In Sec. 
II, we introduce the model Hamiltonian for broad res- 
onances and study the formation of dimers in a static 
Fermi sea. We analyze the effects of Pauli blocking on 
the dispersion of dimers, and illustrate that a kinematic 
effect in the limit of very light minority atoms favors 
finite-momentum dimers near broad resonances. We then 
present the results for quantum mixtures so far studied in 
laboratories. In Sec. Ill, we go beyond the static-Fermi- 
sea approximation to take into account various correc- 
tions due to particle-hole fluctuations near the Fermi sur- 
face of majority fermions. Our calculations indicate that 
the main effect of fluctuating pairs is to reduce the bind- 
ing energy of dimers and also to disfavor the formation 
of dimers with total momenta close to Kkp. In Sec. IV, 
we discuss the limitations of our results, explore the im- 
plications on near-resonance quantum mixtures with a 
finite density of minority atoms and comment on various 
results 0, [TH, [H, |H, obtained in previous studies. 
We summarize the results and discussions in Sec. V. 



II. ANOMALOUS DIMERS: PAULI BLOCKING 
EFFECT 



We start with a single-channel Hamiltonian with a 
short-range interaction Ubf that was introduced previ- 
ously for a study of Fermi-Bose mixtures near broad res- 
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Here /^^(/kj&k) are creation (annihilation) operators 
for Fermi and Bose atoms respectively, and e^ B ^ — 
Jt) k \ — are kinetic energies for fermions (bosons) and fl is 

the volume. Ubf is the strength of short-range interaction 
between fermions and bosons, and is equal to scattering 
lengths 2Trh 2 a/mR after regularization [39( (see also Ap- 
pendix A) . And iriR = mBTUF / {tub +mp) is the reduced 
mass, = ft 2 |k| /2m,R. This single-channel model effec- 
tively describes near-resonance physics provided the res- 
onance is broad enough; i.e., the effective range is much 
smaller than the Fermi wavelength [49j . 
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FIG. 1: (Color online) (a)-(c) Schematics of the Pauli 
blocking effect on dimers with total momentum |fiQ| = 0, 
0.2(m_g /mn)hkF , (ms/m^hkF, respectively. Yellow (light 
gray) areas in the relative momentum k r space are for the 
blocked or occupied states and blue (dark gray) areas stand 
for states available for pairing near the threshold of the two- 
body continuum, (d) Energy dispersion of excited dimers for 
the mass ratio ros/mF = 0.05; shown in the inset is the 
part of dispersion for |Q| > &f. Shaded areas are for the 
two-atom excitation continuum. At l/fc^a = 0.326, the mini- 
mum of dispersion reaches zero energy at a finite momentum 
foQmin = 1.23hkF- Beyond this point, a dimer should appear 
in the ground state when a minority atom nearly resonates 
with a Fermi sea. 

To get insight into pairing with a nonzero total mo- 
mentum ?iQ, we examine the dispersion for a two-body 
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m^m B *Q/k F 

FIG. 2: (Color online) For a dimer with total momentum hQ 
to appear in the two-body excitation spectrum, l/(fcpa) has 
to exceed a minimum value indicated by the curve. Dashed 
line is for free space without a Fermi sea. 



bound state of Fermi-Bose atoms with an arbitrary to- 
tal momentum HQ or a center-of-motion kinetic energy 
£q = h 2 Q 2 /2rriT and wit = mp + tub- The dispersion 
can be obtained by solving the two-body problem on top 
of a "frozen" Fermi sea of majority atoms. The resultant 
self-consistent equation is 
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Here Wb (Q) is the energy of dimers with momentum TiQ, 
and is measured from the Fermi energy ep — H 2 kp/2mF 
of majority atoms. The unit step function 0(- • • ) in the 
sum excludes occupied states below the Fermi surface. 

This equation can also be obtained by studying the 
scattering matrix in the presence of a Fermi sea; the pole 
of the T matrix on the real axis, if it exists, corresponds 
to a bound state with an infinite lifetime (see Appendix 
A). We find that Eqj2]is very similar to the well-known 
Cooper's solution to two attractive electrons on the top 
of a Fermi sea[5(|. Note that the T matrix approach 
usually yields an additional contribution to the Cooper's 
original solution due to the inclusion of hole like config- 
urations in the eigenvalue equation. However, since here 
we are dealing with a situation where only a single mi- 
nority atom resonates with majority ones and there is no 
minority Fermi sea, in our case the hole like configura- 
tions do not contribute to the binding energy. 

The presence of a Fermi surface leads to an anomalous 
dimer excitation spectrum when the scattering length is 
small and negative. Two interesting aspects of the dis- 
persion are worth emphasizing. First, dimers with zero 
momentum appear in the form of excitations for arbitrary 
negative scattering lengths even when they are small in 
magnitude. On the contrary, dimers with a finite mo- 
mentum can exist only when the scattering lengths ex- 
ceed a critical value. In fact, due to the Pauli blocking 
effect, the threshold of two-body continuum for a total 
momentum /iQ is 
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different from the case in vacuum where 
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As a result for Q = 0, one finds that 



the available density of states D2b{E) for two-body scat- 
tering does not vanish when E approaches E t h- When 
< |Q| < ^-kp, one finds that D2b(E) always vanishes 
linearly as a function of energy E — E t h, i.e. 
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The qualitative difference between D2b(E) for the zero- 
momentum channel and for the finite-Q channels is 
schematically illustrated in Figs. 1(a)— (c). Note that 
near the threshold, the two-atom scattering states with 
zero total momentum Q = are represented by the whole 
shell around the Fermi surface while the scattering states 
with finite total momenta only correspond to a small strip 
around the Fermi surface. The relatively low density 
of states near E t h for finite-Q scattering puts a severe 
constraint on the formation of dimers with a finite mo- 
mentum. And for a given negative scattering length, we 
find that dimers only appear in the spectrum up to a 
maximum total momentum. This is reflected in Fig. 1(d) 
in which the dimer dispersion ends up at a finite mo- 
mentum and merges into the two-atom scattering con- 
tinuum. Alternatively, one can conclude that to form a 
finite-momentum molecule, the value of (fcpa) -1 has to 
exceed a minimum value as shown in Fig. 2. 

Second, the dispersion minimum might locate at a 
finite center-of-mass momentum ?iQ. And this is the 
case either when atoms are away from resonances with 
small negative scattering lengths or when Bose atoms are 
very light. In the former case, it is due to the decreas- 
ing threshold E t h as the momentum Q increases from 
zero and therefore the energy of bound states below the 
threshold (see Fig. 1). In the latter situation, to form a 
molecule with zero total momentum HQ = Bose atoms 
need to be promoted at least to right above the Fermi 
surface resulting in a high energy penalty e^, while to 
form a molecule with total momentum near hkp Bose 
atoms need not to be elevated. So energetically it could 
be more favorable to have molecules of Fermi-Bose atoms 
with a nonzero total momentum /iQ in this limit. 

We have examined m e f / , the effective mass of extended 
dimers near Q = 0, as a function of scattering length 
a and the mass ratio ms/rnp. At any small negative 
scattering length [— (kpa)^ 1 3> 1] or far away from reso- 
nances, 
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and it is indeed always negative. At scattering lengths 
where creation of a dimer with zero momentum costs 
no energy or Wb(Q = 0) = 0, we find that 



1 



l eff 



1 

rriT 



Amp 
3m n' 



m R 
mp 



(6) 



4 




FIG. 3: (Color online) (a) Critical values of a at which the 
dimer excitation energy Wb(Q) becomes zero at Q — Qmin , 
the minimum of the excitation spectrum; a dimer starts 
emerging in the ground state for a minority atom resonat- 
ing with a Fermi sea beyond this line. The dashed line is 
obtained by setting Wb{0) = or for the scattering lengths 
at which the zero-momentum dimer energy vanishes, (c) is 
for larger mass ratios, (b) Qmin, the momentum at which 
the dimer dispersion touches zero for the first time when ap- 
proaching the resonance from the side of negative scattering 
lengths, as a function of mass ratio ms/mi?. These plots are 
obtained assuming the Fermi surface is static and there are 
no particle-hole fluctuations. 



and the dimensionless function 
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Note that m e f f in Eq|5] is an indicator of relevance or 
irrelevance of zero-momentum dimcrs when dimers start 
appearing in the ground state near resonance. And as 
far as ms/rap > 0.7 and the energy penalty €p is not 
too heavy, m e // is positive, although it can be much big- 
ger than the bare mass rar(= nif as a result of 
dressing in the Fermi sea. In this limit, the dimer spec- 
trum indeed crosses zero first at Q = and a zero mo- 
mentum dimer is expected to appear in the ground state. 
However, when nig < 0.7mp, the effective mass becomes 
negative, implying the relevance of dimers with finite mo- 
menta in the reconstructed ground state. Indeed we find 
that for small mass ratios, the minimum of the dispersion 
spectrum crosses zero first at a finite momentum hQ m i n 
when approaching resonances from the side of negative 
scattering lengths. Going further beyond this point, one 
would expect that dimers with finite momenta hQ m i n 
thus appear in the ground state for minority atoms res- 
onating with a Fermi sea. Details are shown in Fig. 3. In 
Figs. 4-5, we also present the results for dimers near a 



few broad interspecies resonances that have been studied 
in recent experiments. 
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FIG. 4: (Color online) Energy dispersion of dimers near an 
interspecies resonance of 6 Li- 40 K atoms. ms/m_F = 0.15 for 
a 6 Li atom nearly resonating with a Fermi sea of heavy 40 K 
atoms. Shaded areas are again for two-atom excitation contin- 
uum. At 1/kFd = 0.393, the minimum of dispersion reaches 
zero energy at a finite wave vector Qmin = 1.05&;f. 



III. EFFECT OF FLUCTUATING 
PARTICLE-HOLE PAIRS NEAR THE FERMI 
SURFACE 

Besides the Pauli blocking effects, the energetics of 
bound states can be further modified by the dynamics of 
Fermi seas, or virtual particle- hole pairs near the Fermi 
surfaces, so the dimers are further dressed in fluctuating 
particle-hole pairs. One of the main effects of these fluc- 
tuating pairs is to renormalize the two-body interactions 
near the Fermi surface and therefore to modify the bind- 
ing energetics of dimers as well as the dispersion. The 
fluctuating particle-hole pairs also lead to a decay of the 
dimer into either another lower energy dimer or two un- 
bound atoms in the continuum and we will not discuss 
the decay in this paper [5lT ] . 

The main effect on two-body scattering and therefore 
the dimer binding energy is illustrated in FiglHl both 
schematically and diagrammatically. Consider a general 
situation where the density for minority atoms or Bose 
atoms is finite but much smaller than the density of ma- 
jority Fermi atoms [52j. Three classes of particle-hole 
pair fluctuations contribute: 

(A) A majority Fermi atom is first excited from be- 
low the Fermi surface to above by the incoming colliding 
minority Bose atom; the Fermi hole is later filled by the 
incoming Fermi atom that simultaneously scatters the 
minority Bose atom to its final state. 

(B) A minority Bose atom is excited from the conden- 
sate by the incoming colliding majority Fermi atom; the 
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FIG. 5: (Color online) Energy dispersion of dimers for dif- 
ferent mass ratios, (a) ms/mF = 6.67 for a single 40 K 
atom resonating with a Fermi sea of light 6 Li atoms; (b) 
ms/mF = 2.175 or for an interspecies resonance between 40 K 
and 87 Rb atoms; (c) ms/mF = 23/6 or 6 Li- 23 Na resonance; 
(d) ms/mF = 14.83 or 6 Li- 87 Rb resonance. In (b)-(d), mi- 
nority atoms are Rb, Na, and Rb atoms, respectively. Shaded 
areas are again for two-atom excitation continuum. In (a)- 
(d), the minimum of dispersion first reaches zero energy at 
zero momentum. 



Bose hole is later filled by the incoming Bose atom that 
simultaneously scatters the Fermi atom to its final state. 

(C) A majority Fermi and a condensed Bose atom are 
excited with holes left behind filled in later by the incom- 
ing majority Fermi and minority Bose atoms. 

A- and C-class processes involve exchange of virtual 
Fermi particle-hole pairs; i.e., a second majority Fermi 
atom has to be ejected out of the Fermi sea; they both 
have an additional negative sign with respect to the pro- 
cesses in class B (see Fig. [6]) or the direct processes with- 
out virtual particle-hole pairs. One can also carry out a 
parallel analysis on a Fermi-Fermi mixture by replacing 
a condensate with a Fermi sea of minority atoms. In that 
case, the processes in the A-class or B-class will have a 
different sign with respect to the processes in the C-class 
as well as the direct scattering. And in the limit of a sin- 
gle minority atom that we are going to focus on or when 
the minority atom density vanishes, only the processes in 
the A-class contribute. 

To understand this particular energetic effect of fluc- 
tuating particle-hole pairs, we study the T matrix in- 
cluding vertex corrections due to these fluctuating pairs. 
The technical details of this part are summarized in Ap- 
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FIG. 6: (Color online) An illustration of induced scattering 
by virtual particle-hole pairs near a Fermi surface (a) , a con- 
densate (b), or both (c). (d) is for direct scattering with- 
out virtual particle-hole pairs. In (a)-(d), a minority atom 
[blue(dark gray) dots] and majority atoms [ yellow(light gray) 
ones] have an interspecies resonance with zero total momen- 
tum; the light blue (light gray) area represents states occupied 
by fermions. Also shown is a diagrammatic representation of 
induced scattering with the time order illustrated explicitly. 
Solid line is for Fermi atoms, dashed line is for minority Bose 
atoms with nonzero momentum, and wavy lines are for Bose 
atoms in a condensate. In (e) and (f), we also illustrate two 
distinct virtual processes in (a)-(d): (a) and (c) involve, in 
addition to the incoming Fermi atom, a hole below the Fermi 
surface [as shown in (e)] while (b) or (d) does not [see (f)]. 
Because of the Fermi exchange statistics, the amplitude of 
the (e) process has an additional negative sign with respect 
to that of the direct process (d) as well as that of the (f) pro- 
cess. In the single-minority-atom limit, only A-class and the 
direct processes contribute; the dashed line can also represent 
a minority Fermi atom in this limit. 



pendix B. Scattering with the Fermi sea effectively takes 
momenta away from the incoming minority atom of the 
dimer to create virtual particle-hole pairs. This virtual 
process crucially depends on the initial and final relative 
momentum between the two atoms in the dimer. This 
can been illustrated via a momentum flow chart in Eq. 
(|5|) for two incoming atoms with the total momentum 
fi,Q and relative momentum Kk. scattered into two outgo- 
ing atoms with the same total momentum but different 
relative momentum fik' through colliding with the Fermi 
surface. 
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where \F.S. > is introduced to represent the Fermi sea of 
majority fermions; the virtual state energy explicitly de- 
pends on Q, k, k'. So with the particle-hole fluctuations, 
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the T matrix not only depends on the total momentum 
S-Q but also the relative momentum hk and Hk' . Effec- 
tively, the fluctuating pairs mediate an interaction with 
a finite range. 

One of the ways to estimate these effects is to introduce 
an effective scattering length (see Appendix B for de- 
tails) . The effective scattering length a can be expressed 
as 

a' 1 = a' 1 - k F R + 0(k 2 F a)..., (9) 

Here a differs from the free space scattering length a due 
to the particle- hole pair fluctuations. And kpR repre- 
sents the lowest order vertex correction induced by fluc- 
tuating particle-hole pairs near a Fermi surface (or in 
condensates if a finite density of minority Bose atoms are 
present, or both) as shown in Fig. [6] (see also Appendix 
B). This correction is analogous to the vertex corrections 
discussed for zero-momentum pairing in fermion super- 
fluids by Gor'kov and Melik-Barkudarov (53|. We shall 
denote this correction as the GMB vertex corrections (see 
also [Hj]) in the following discussion. Here we have car- 
ried out a detailed diagrammatic analysis on the GMB 
correction in finite-Q channels. We obtain the momen- 
tum dependence of induced scatterings and summarize 
the main results below. As a remark, we find that in the 
limit of heavy minority atoms (large mass ratio mj /mp), 
the R function for two-atom scattering with zero total 
momentum is 

R^Uog^. (10) 

7T mp 

We use natural logarithmic functions with base e here 
and throughout the paper. This logarithmic divergence 
in the large-mass-ratio limit can be attributed to the 
very heavy dressing of a zero-momentum dimer in vir- 
tual particle-hole pairs near the Fermi surface. One can 
also demonstrate that this has the same origin as the 
Anderson's infrared catastrophe in a quantum impurity 
problem [54j. 

We have numerically found that for scatterings near 
the threshold of two-particle continuum, the exchange 
processes in the A-class described above always induce an 
effective repulsive interaction. So the fluctuating particle- 
hole pairs in this case effectively screen or reduce the 
interspecies attractive interactions in all channels of mo- 
menta HQ (see Fig(7|) . Furthermore, the magnitude of the 
corrections to the effective two-atom scattering strongly 
depends on the total momentum HQ of the scattering 
atoms. This can be demonstrated by explicitly examin- 
ing two cases: HQ — and HQ = hkp in the limit of light 
minority atoms, i.e., when the energy of virtual states is 
dominated by the minority atom. For a pair of minor- 
ity and majority atoms in the HQ = channel, or with 
momentum (Hk, — Hk), it requires that the minority atom 
momentum be close to hkp because of the Pauli blocking 
effect. Consider such an incoming minority atom with 
momentum Hk colliding with the Fermi surface popping 
up a pair of particle-hole excitation. The intermediate 



virtual state with total energy e v , and total momentum 
Hk consists of an outgoing minority atom and a particle- 
hole pair of majority atoms. This energy e v can be ei- 
ther larger or smaller than the energy of the incoming 
minority atom ■ The contribution to the effective in- 
teractions from a virtual state of energy e v is inversely 
proportional to — e v and so the contributions from dif- 
ferent virtual states would have different signs depending 
on whether is larger or smaller than e v , thus leading 
to a destructive interference between different configura- 
tions. Now let us turn to the scattering channel with the 
total momentum near fi|Q| = hkp, i.e. a minority atom 
with hit = interacting with a majority atom near the 
Fermi surface with Fermi momentum hkp. For a minority 
atom with = 0, all intermediate virtual states would 
have energies larger than the initial one and all should 
contribute to the effective interaction with the same sign, 
leading to a constructive interference. This results in a 
maximum value of the correction to two-body interac- 
tions; i.e., the correction reaches a peak value when hk 
is near zero or fi|Q| is near hkp. 

A more elaborated examination along this line can be 
carried out in Appendix B. There we show that for all vir- 
tual states involving a hole like excitation of momentum 
HI, the most dominating contribution to the overall vertex 
correction is always from the state with HI ~ HkpQ/\Q\ 
so that the energy difference between the virtual state 
and the initial state is minimum. We can further intro- 
duce an effective group velocity for virtual states involv- 
ing a hole like excitation with momentum as HI, de v (1) /d\. 
For the scattering near the threshold of the two-body 
continuum, this velocity turns out to be proportional to 
mpi/rriBQ — kp. In the limit of light minority atoms, 
i.e., the reduced mass ma approaches mg, the dominat- 
ing virtual states appear softest near HQ = hkp, consis- 
tent with the above qualitative argument based on in- 
terference effects. In general, formation of bound states 
near HQ = {tub I rnp)hkp is disfavored by the fluctuating 
particle-hole pairs. For this reason, the vertex correction 
to two-body interactions due to the creation of particle- 
hole pairs is relatively small near small and large total 
momentum, and peaks near HQ — Hkp in the limit of 
light minority atoms. This is confirmed in our numerical 
evaluations of R. So fluctuating particle-hole pairs fa- 
vor the formation of dressed bound states with HQ away 
from hkp, over the dimers with HQ ~ Hkp. Note that 
although we arrive at this conclusion by considering only 
the lowest order corrections, we conjecture that this is 
generically true as far as the effective group velocity for 
the hole like excitations has a minimum near HQ = Hkp. 
The position of the peak in FigfT] in our opinion corre- 
lates with the threshold minimum where the screening 
effect is usually the strongest. 

The problem of two-atom scattering near HQ = hkp 
is closely related to the spin polaron in cold gases 
that has been discussed quite extensively in the recent 
literature [3ll - [38| . In our calculation, we also see a transi- 
tion from a polaron to dimer when approaching the res- 



7 



onance from the negative-scattering-length side. We re- 
fer the readers to those publications for more elaborated 
discussions on the spin polaron. However, we should 
also emphasize that although the physics of the spin po- 
laron itself is an interesting subject, to understand the 
nature of the ground state of a single minority atom; 
or the general pairing correlations in quantum mixtures 
near resonances, it is necessary to study the whole dimer 
spectrum. This is because as we have seen before, the 
minimum of the dispersion can occur at any momentum 
depending on the mass ratio ^s.. The minimum loca- 
tion hQmin generally differs from HQ = hkp or HQ = 0. 
Only for a certain range of mass ratios, scatterings near 
HQ = or HQ = kp set the overall instability of a un- 
paired minority atom. 
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(b) K40-Rb87 
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FIG. 7: (Color online) (a) R, the vertex correction as a func- 
tion of total momentum HQ for mass ratio ma/mf = 0.135. 
Here the parameter 7? is defined as 1/ka — — R. Dif- 

ferent plots are obtained with different choices of energy E; 
the solid line is obtained using a self-consistent method where 
E = Wb(Q)- (See more discussions in Appendix B.) (b) Dis- 
persion of dimers with GMB and without MF particle-hole 
fluctuations or formally Gorkov-Melik-Barkudarov (GMB) 
vertex corrections for ms/mF = 0.135 and (fepa) 1 = 0.62. 
(c) R for the Q — scattering channel as a function of mass 
ratio ms/niF. 



IV. IMPLICATIONS 

The dispersion obtained for an individual minority 
atom having interspecies resonance with majority Fermi 
atoms that form a Fermi sea can be applied to identify the 
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FIG. 8: (Color online) The vertex corrections to the dimer 
spectrum (denoted as GMB). Thin solid lines are the refer- 
ence spectrum without the vertex correction, or a mean field 
(MF) result, (a) tub/thf = 6.67 for a single 40 K atom res- 
onating with a Fermi sea of light 6 Li atoms; (b)m,B /niF = 
2.175 or for an interspecies resonance of 40 K- 87 Rb atoms; (c) 
tub/tiif = 23/6 or 6 Li- 23 Na resonance; (&)m,B /rriF = 14.83 
or 6 Li- 87 Rb resonance. In (b)-(d), minority atoms are Rb, 
Na, and Rb atoms respectively. Shaded areas are again for 
two-atom excitation continuum. 



instability point [551] for an extremely imbalanced quan- 
tum Fermi-Bose or Fermi-Fermi mixture beyond which 
the weakly interacting mixtures become unstable. Be- 
low we assume the minority species can be either Fermi 
or Bose atoms with densities much lower than that of 
majority atoms. 

The instability point corresponds to the scattering 
length at which the dispersion of anomalous dimers 
touches zero or the energy cost of creating a dimer be- 
comes zero. Since the minimum in the dimer dispersion 
can be either at Q = or Q ^ 0, it is implied that the 
instability might be driven by formation of dimers with 
either Q = or a finite Q. 

We summarize the results of the dispersion for differ- 
ent mixtures in Fig. 8, and the instability for different 
mass ratios in Fig. 9 Since the fluctuating particle-hole 
pairs near the Fermi surface reduce the binding energy of 
dimers, the formation of a dimer in the ground state for a 
minority atom and a Fermi sea of majority atoms occurs 
at a higher critical value of \jkpa when crossing the res- 
onance from the side of negative scattering lengths. The 
second effect of fluctuations here appears to destabilize 
dimers with finite momenta near kp , because of the peak 
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structure in the vertex corrections. This is qualitatively 
consistent with the results in Fig. 9(b): The vertex cor- 
rections suppress the area of the region where Q m in is 
finite. 

For the case of equal-mass interspecies resonance, the 
mean field calculation with a static Fermi surface yields 
a critical value, (kpa)^ 1 — 0.24; fluctuations of particle- 
hole pairs correct the mean field instability point lead- 
ing to a modified critical value (kpa)^ 1 — 0.34. Com- 
pared to (fcpa) -1 = 0.88, the critical value obtained via 
the diagrammatic quantum Monte Carlo [34j that has so 
far been done only for the equal-mass case, our approx- 
imation yields a qualitative correct account of the main 
effects of fluctuations but quantitatively underestimates 
the effects of fluctuations. And in our diagrammatic ap- 
proach, for mass ratios (ms/ra^) less than 0.135 a dimer 
with finite Q appears in the ground state for a minority 
atom resonating with a Fermi sea of majority atoms at a 
critical value of scattering length a. This is close to 0.15, 
the numerical result obtained in a truncated particle-hole 
pair subspace [liij . However for this critical mass, the 
critical scattering length at which a dimer starts forming 
in the ground state is predicted as (kFa)~ l — 0.62 in our 
approach while in Ref. 18], (kpa)^ 1 is close to a larger 
value of 1.7. 



(a) (b) 




0.2 0.4 0.6 0.2 0.4 0.6 

m B /m F m B^ m F 

FIG. 9: (Color online) The estimated effect of particle-hole 
fluctuations(curve labeled as GMB) on the dimer formation; 
the result obtained in the static Fermi surface approximation 
is denoted as MF. (a) Critical line beyond which a dimer starts 
to appear in the ground state for a minority atom resonating 
with a Fermi sea of atoms; the dashed lines are obtained when 
only taking into account the zero-momentum dimers, i.e. set 
by scattering lengths at which Ws(Q = 0) vanishes, (b) The 
momentum of a dimer that appears in the ground state right 
above the critical line. 

It remains to be understood what happens beyond the 
instability line identified in FigO and what is the na- 
ture of the transitions if there are any [56j . For Fermi- 
Bose mixtures, an earlier attempt was made to under- 



stand dimer correlations beyond a critical line [29[ that 
was obtained when the scatterings in finite- Q channels 
were not included. However, the appearance of finite- 
density fermion dimers in the ground state beyond the 
critical line was not properly taken into account in anal- 
yses. This complication and more importantly potential 
dimer-dimer and atom-dimer interactions about which 
we know very little make a thorough study in this limit 
extremely challenging. The second issue is how dimers 
in a finite- |Q| channel compete against other structures 
with higher order correlations. For instance, it was also 
proposed that for two-dimensional Fermi- Fermi mixtures 
with extremely small mass ratios less than 0.01 [15j, a 
crystal structure can further develop at exponentially low 
temperatures. This question perhaps can be best ad- 
dressed by numerically taking into account various fluctu- 
ations. On the other hand, we remark that LOFF states 
have been suggested in Fermi gases near resonances. For 
Fermi-Fermi mixtures with equal masses, Bulgac et al. 
suggested that the ground state can be a pairing state 
that breaks translational symmetries [l4j . Recently, it 
was also argued that LOFF states appear on the molec- 
ular side of resonances in almost fully polarized Fermi 
gases but only for nonequal masses; in addition, a trimer 
phase was also proposed [l8|. It remains to be clarified 
the connections between different results obtained in Ref. 
[H, Ref. [3 and Ref. [H[. In Ref. [11], the authors 
first investigated the anomalous dispersion of molecules 
in Fermi-Bose mixtures and estimated the effect of fluctu- 
ating particle-hole pairs on the formation of dimers. The 
authors then proposed a first-order phase transition be- 
tween the noninteracting ground state and a fully paired 
state analogous to the BCS one. Dukelsky et al. [57| 
correctly pointed out that a string factor was missing in 
the energetic analysis of the first-order phase transition. 
The correct form of the energy function should be 

e= E(ki 2 + i^i 2 )^ + Eh 2 ^ 

k k 

, Ubf * * 

+— 2^ u k'«k' u k«k x 

kk' 

n (K-i 2 -K"i 2 -k»i 2 )- in) 

k"e(k,k') 

Here the string product is carried over states between k 
and k' in the trial wave function given in Eq.(5) in Ref. 
(391 ] - In the previous calculations, the string product was 
indeed overlooked and its effect was not properly taken 
into account. The string factor modifies the interaction 
energy via inducing a rapidly alternating sign when the 
momentum varies, and more importantly unless the dis- 
tribution function is strictly the Fermi-Dirac function, 
the string product should be equal to zero in the ther- 
modynamic limit where there are an infinite number of 
states between the two open ends of the string. For this 
reason, the energy of the proposed pairing states near 
the first-order phase transition was erroneously underes- 
timated and within this mean field approximation with 
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the string product properly taken into account, one is 
no longer able to show that there is a first-order phase 
transition to the proposed pairing state as argued in the 
previous Letter [39| . In Ref . [H?} , the Richardson solution 
has also been further utilized to argue the nonexistence 
of collective states below the energy of zero-momentum 
molecules. It was pointed out that Fermi-Bose molecules 
are noninteracting in the reduced pairing model employed 
there; in that reduced Hamiltonian that differs from the 
full scattering Hamiltonian [see Eq.(l)] introduced here, 
scatterings in all finite-Q sectors are completely sup- 
pressed. So the result of the nonexistence of lower col- 
lective modes, although it is a feature of the reduced 
Hamiltonian, might not be a property of cold atoms 
near Feshbach resonances. In fact, one can show that 
the second-order commutator-anticommutator for com- 
posite molecules does not vanish when the full Hamilto- 
nian in Eq.Q) rather than the reduced pairing model is 
considered [58[. Two-body scattering with a finite to- 
tal momentum HQ that was neglected in the reduced 
pairing Hamiltonian induces significant dimer-atom or 
dimer-dimer scatterings as well as scatterings between 
a dimer and a Fermi sea. In Fermi-Fermi mixtures, the 
finite-Q scattering also plays a critical role in dimer-atom 
or dimer-dimer scattering near resonances. The dimer 
or molecule dynamics near resonance can be understood 
when these scatterings are properly taken into account. 
For this reason, the issue of low-lying collective modes of 
the molecules remains an open question and can be ad- 
dressed when the physics beyond the Richardson-solution 
is considered. 

In summary, to understand the physics beyond the in- 
stability line, it is essential to include the two-body scat- 
tering with a finite total momentum fi,Q when locating 
the critical scattering length for the molecule dispersion 
can be quite anomalous and its mass near Q = can be 
negative either when far away from resonance or when 
the minority atom is very light. In the latter case, one 
can show that Fermi-Bose molecules with finite momenta 
should first appear in the ground state because Q = 
molecules are energetically unfavorable; the instability 
of the noninteracting ground state is therefore driven 
by scatterings in finite-Q sectors rather than two-body 
scattering processes in the reduced pairing Hamiltonian, 
which excludes two- atom scattering with finite Q. 

V. CONCLUSION 

In conclusion, we have investigated dressed dimers in 
quantum mixtures and carried out thorough studies of 
energies of dimers in finite-Q channels. Our calculation 



suggests that in the limit of light bosons, dimers with a 
finite momentum HQ are most relevant for the ground 
state for a minority atom resonating with a Fermi sea 
of majority atoms. There are two open issues that need 
to be understood in the future, most likely beyond the 
framework of the approach here. One is the role of three- 
body correlations in quantum mixtures. Following a gen- 
eral argument by Efimov, when ms/mF mass ratios are 
smaller than 0.145, trimers involvin g tw o heavy majority 
Fermi atoms should form in vacuum [59j . However, three- 
body dimer-atom resonances induced by trimers gener- 
ally are much narrower than two-body resonances and 
contributions from Efimov physics to the two-body effec- 
tive interactions can be smaller than the universal correc- 
tions that originate from collective excitations [53| which 
we have taken into account in this article. This was pre- 
viously noticed in the studies of Fermi-Fermi mixtures 
[l6j]. In this limit, one should expect that trimers have 
relatively weak effects on the physics in two-body chan- 
nels and that dominating effects in the dimmer channel 
discussed above can survive these higher order correc- 
tions. 

A related and perhaps more important question is un- 
der which conditions three-body correlations are domi- 
nating and an additional trimer channel has to be added 
to the discussion on many-body physics. How trimers 
mediate additional strong scattering between dimers or 
dimer and atoms and lead to nontrivial higher order 
many-body correlations represents one of the most chal- 
lenging issues we face now and remains to be studied 
in the future. The answer to this question should de- 
pend on short-distance behaviors of Efimov potentials 
and would therefore be nonuniversal. Our results about 
dimers should be relevant whenever these nonuniversal 
higher order correlations are insignificant, and the frac- 
tion of atoms forming trimers is small. 
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FIG. 10: Feynman diagrams for the T matrix calculation, thick (thin) solid lines are for majority (minority) atoms, dashed 
lines are for the bare interaction Ubf, and T(- • ■ ) is the irreducible two-body interaction vertex, (a) The ladder summation and 
diagrammatic representation of the Bethe-Salpeter equation in Eq. (|A4I) . (b) T with GMB corrections. Here the interaction 
lines within the Fc(l, 2) are renormalized interactions, which can be expanded in series of ladder diagrams using bare interaction 
U bf . 

Appendix A: T matrix in the static-Fermi-sea approximation and the molecule's effective mass 



To determine the molecular excitation energy in a Fermi sea of majority atoms, we look at the T matrix for the 
minority-majority scattering. Diagrammatically it is the skeleton diagram (without the external propagators) of the 
following time-ordered two-particle propagator: 



(F.S.\Tc(t, k 2 )b(t, Q - k a )&t(0, Q - kx) c t(0, k^F.S.) 



(Al) 



where \F.S.) stands for the Fermi sea of majority atoms. An isolated pole of T in the energy space represents a dimer 
excitation. The pole could be either on the real axis or in the lower complex half plane; in the latter case the dimer 
has a finite lifetime. Note that the above time-ordered propagator is zero when t < 0, because one cannot create a 
hole excitation of minority atoms. 

In the following we use the ladder approximation to calculate T as shown in Fig llOl The rung of each ladder is 
the r function or the irreducible two-body interaction vertex, and the side rails are Green's functions of majority 
and minority atoms. One can write down the Bethe-Salpeter equation to include all the ladder-like diagrams as the 
following: 



T(l,2) = T(wi,ki,£-wi,Q-ki;w2,k2,E-W2,Q-k 2 ) 
ff(3) = iG F (uj 3 ,k 3 )GB{E-Lo 3 ,Q-k 3 ) 

T(l,2) = r(l,2) + J]r(l,3)fl(3)T(3,2) 



(A2) 
(A3) 
(A4) 



where we have introduced a shorthand notation (1,2) to denote the incoming and outing states. Here Gf (Gb) 
are bare Green's functions for majority(minority) atoms defined as, G^ 1 (oj,k) = uj — + i0 + and G i ^ 1 (aj,k) = 
uj — (e£ — e f")(l — *0 + ). In the leading order approximation, we have r(l, 2) = Ubf, and one can further integrate over 
w 3 . The result is 
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T(l,2) 
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Wf 



n 4^ 



l-np(k) 



E-i0+ 



(A5) 
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The energy of a dimer is determined from the condition that T becomes divergent when E approaches Wb(Q)- 



2Trh 2 a 



n 



1 -n F (k) 



Q-k 



W B (Q) 



(A6) 



The last term in Eq. (|A6j) was introduced to regularize the ultraviolet divergence. 

The integral in Eq. (IA6|) can be evaluated analytically for an arbitrary value of Wb{Q)- And the result for the 
bound-state solution when Wb(Q) < E t h(Q) is 
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1 



2b + c 
2/n2 



+ \Jb 2 - c 
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(A7) 



2ms 



~ e£ ^ 2m B 



Wi 



11 



From the above expression one can find the effective mass of the molecule m e // defined as Wb(Q) — Wb(0) — 2m > 

1 1 2ni R 1 1 

m e ff m B + m F 3m| 1 - y 2 1 + ±^jt log ±±| 

» - ^'W = °» + ^ (A8) 

The effective mass at the instability point given in Eq. © is obtained by setting Wb{Q = 0) = 0orj/ = x= y^J- 
One can also find out that the effective mass is always negative in the limit where ^— - — » — oo by studying the 
asymptotics when y — > 1 — 2 exp f ^|^J : 

1 TOfl / 7T 



m e // 6m^ \ kpa 



2 ^p ( -— ) (A9) 



Appendix B: Effective scattering length a From T matrix Method 

In this section we use the T matrix method to study the effects of fluctuating particle-hole pairs near Fermi surface. 
Our analysis shows that in the limit of small kpa, the effects due to particle-hole fluctuations can be captured by 
introducing an effective scattering length a as 5 _1 = a -1 — kpR + kpO(kpa), where the function R represents the 
lowest order vertex correction. Similar quantities such as effective interactions have been introduced before 
to address effects of particle-hole fluctuations in zero-momentum pairing physics. The T matrix method used here 
provides a unified description for dimers of different momenta and energies, and can be further extended to near 
resonances (5lT |. In the following we first show how to obtain Eq. (|B7[) for dimer's energy by examining the pole 
structure of the T matrix. Then we show how to obtain Eq. (|B15P where an effective scattering length emerges. 

We start from Fig fTCT b). where we include the leading order term Fc7(l,2) as 

r(l,2) ~ U bf +T c (l,2) (Bl) 



r c (l,2) ~ / |jMM)GB(w-w 1 +i5- Wa ,Q-k 1 -k a 



1) 



2_7r^a\ l_y^ n F {\) 

m R J E _ L0l _ Lj2 _ e F +e F _ e B +1 _ ki _ ka+i 0+ 

Here Tc(l, 2) corresponds to exactly the type-A virtual process discussed in the main text. Note that in Tc(l, 2) we 
replace the renormalized interaction vertex (which is a sum of diagrams with bare interactions) for particles or holes 
near the Fermi-surface by two-body scattering amplitude in vacuum, i.e., 2 ^ a - This is justified in the limit of small 
kpa; such a replacement will not change function R, but only leads to inaccuracies in higher order terms. With Eq. 
(|A4[) one can write down the corresponding Bcthc-Salpctcr equation: 

T(l,2) = [[/ b/ + r c (l,2)]+^[[/ t/ + r c (l,3p(3)T(3,2) (B3) 

3 

The above equation defines an infinite series implicitly in terms of Ubf and Tc- It is not hard to regroup terms in the 
series according to the powers of Ubf and sum them up in the following way as 



T(l,2) 



<Mi)<M2) 



>M 



U b -/ - E 3 H (3) + E 34 H (3)r c (3, A)H (4) + E 345 ^ (3)r c (3, 4)H (4)r c (4, 5)H(5) + ■■■} 

r c (i, 2) + r c (i, 3)ff(3)r c (3, 2) + ^ r c (i, 3)ff(3)r c (3, 4) J ff(4)r c (4, 2) + • • • I (B4) 

3 34 J 

(1) = 1 + ^r c (l,3)i?(3) + ^r c (l,3)fT(3)r c (3, 4) + J2 T c (l,3)ff(3)r c (3, 4)fT(4)r c (4,5) + ••• (B5) 
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We can assume that the series in the curly brackets and 4>m are convergent, based on the estimation that each term 
is at most of order (kpa) compared the preceding term. Then the pole in the T matrix is completely determined by 
the denominator, and the dimer's energy E should satisfy the following equation: 
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34 
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(B6) 
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(B7) 



Here in Eq. (|B7[) we have only kept leading order terms, and we have further integrated over frequencies W3 and 0J4,. 

Compared to Eq. (IA5|) for the static Fermi sea, the above equation is very similar except for an additional triple- 
momentum integral due to particle-hole fluctuations. We start to analyze this triple integral by examining the range 



of various denominators. We first notice that el 



e Q+l-k-p 



E represents the energy difference 



between the dimer and virtual processes, and hence can be either positive or negative depending on the values of 
k, p, 1 (one exception is at E = 0, the denominator is always positive for all possible k,p,l). On the contrary, the 



other two denominators, i.e., e 



Q 



_ k - e£ - E and e F 



eq_ p — e F — E, are always positive for all possible values 



of k, p, which is guaranteed by the definition of bounded dimers; furthermore in the limit kpa — > — , these two 
denominators can be very close to zero as the binding energy is small. These unique features suggest that the value 
of the triple integral mainly depends on the value of X)|i|</c F ( e k + e p — € \ ~ e F + e Q+i-k-p — ^ ~ *0 + ) _1 when k and 
p are near the two-body continuum threshold; i.e., |k| = |p| = kp (so-called "back-to-back scattering") for Q = 0, 
and k ~ p ~ kpQ/Q (so-called "forward scattering") for finite Q. Based on these considerations, one can introduce 
a step function approximation to simplify the dependence over k and p, 



-r- 



n F (l) 



" ! e k ^ e p e l e F + e Q+l-k-p 



E-i0+ 



f c, k F < |k| , |p| < A c 
0, otherwise 



(B8) 



Here Tc is a constant independent of p and 1. It is natural to choose Tc to be the value for "back-to-back scatterings" 
or "forward scatterings" as 



T c (Q = 0,E) 



dO„ a dtt nD 1 



n F (l) 



4tt 4tt Q ^ e F - e F + e B , . -E- i& 
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n F {\) 



e F + e B 



Q • 1 2fcj -if 



-E-i0+ 



(B9) 
(BIO) 



where n p and n q are unit vectors. The effective cutoff Ac specifies the k or I dependence of the integral in Eq. (1B8[) . 
and is estimated to be a few kps in our case. In the following we will see that the specific value of Ac does not enter 
into the leading order correction. The triple-momentum integral thus can be simplified as 



= (-fa) 
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2irh 2 a 
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e(|k|-A c ) 
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(Bll) 

(B12) 

(B13) 
(B14) 



Here Q(x) is the unit step function. In Eq. (|B12|) we approximate E to be the mean field solution so that the first 
sum is — 3"^ according to Eq. (|A6|) . The second term in Eq. (|B12|) is estimated of m ^ F , if (A c — k 2 F ) m^ 1 is large 
compared to the binding energy. 
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With Eq. (|B12p . Eq. (|B13|) and Eq. (|B7|) . one arrives at the main result of this section 

mji 1 v / 1 — np(k) 1 



l - = -y(- 



2Trh 2 a n^\e% + eg_ fe - e£ - E - i0+ 



T C (Q,E) (B15) 



where the effective scattering length a and the function R are 

1 1 2nh 2 



kpa kpa mRkp 



Ref c(Q, E) (B16) 



R = ^^Ref C (Q,E) (B17) 
m R k F 

Several remarks about the effective lengths. First Tc is usually a complex number, with its real value leading to 
dimers' energy shift, and imaginary part relating to dimers' decay or finite lifetime. Here we only include the real 
part to obtain the energetics. Second, effective scattering lengths usually depend on energy and momentum of the 
pair, and we have plotted the R function for different Q and E in FigJT] In the kpa — > 0~ limit when dimers are 
shallow, one can neglect the energy dependence and simply set R — R(Q,E = E t h(Q)). To extrapolate our result to 
near resonances, we obtain Wb(Q) by solving Eq. (IB15|) self-consistently, and the results are shown in Fig(7{b) and 
Figl 

In the limit of large kpa, the above ansatz in Eq. (|B16[) remains correct, although the estimate of T should include 
many other higher order vertex corrections which are not included in the irreducible diagrams shown in Fig. 10 and 
is therefore more involved. Nevertheless, we expect that the diagrams in Fig. 10 capture the most relevant qualitative 
aspect of particle-hole fluctuations. We therefore extrapolate this result to the unitary limit and apply it to atoms 
near interspecies resonance in Sec. III. Note also that for the dimer energetics in zero- or small-Q channels, the vertex 
correction appears in the lowest order and is the most dominating effect of particle-hole fluctuations. In this limit, 
the self-energy effect (i.e., the effective mass, the residue of the Green's function) of the minority particle appears in a 
higher order in terms of kpa. The renormalized chemical potential does not play an important part in the discussion 
of the dimer; it gives an overall shift of the two-body continuum as well as the dimer's energy. We plan to study 
high-order effects in the future. 
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